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The computational prediction of ligand entry and egress 
paths in proteins has become an emerging topic in com- 
putational biology due to the potential for estimating 
kinetic properties of drug binding. These properties are 
related to important pharmacological quantities such as 
the k on and k off rate of drugs [1,2] .We have investigated 
the influence of protein flexibility on tunnel prediction 
using geometric methods by comparing tunnels identified 
in static structures with those found in structural ensem- 
bles of three CYP isozymes. We found drastic differences 
between tunnels predicted in the crystal structures as 
opposed to those predicted in the ensembles [3]. Further- 
more, we found significant differences between tunnels 
identified in the apo versus the holo protein ensembles [3]. 

While geometric prediction provides a good starting 
point for tunnel prediction, in order to estimate kinetic 
properties, more detailed investigations of the ligand bind- 
ing process are required. We have developed a tunnel pre- 
diction methodology, IterTunnel, which predicts tunnels 
in proteins and estimates the free energy of ligand unbind- 
ing using a combination of geometric tunnel prediction 
with steered molecular dynamics and umbrella sampling 
[4]. Applying this new method to cytochrome P450 2B6 
(CYP2B6), we demonstrate that the ligand itself plays an 
important role in reshaping tunnels as it traverses through 
a protein. This process results in the exposure of new tun- 
nels and the closure of pre-existing tunnels as the ligand 
migrates from the active site. We found that many of the 
tunnels that are exposed due to ligand-induced conforma- 
tional changes are amongst the most energetically favor- 
able tunnels for ligand egress in CYP2B6 [4] . 
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